Interpretation of Quasielastic Scattering Spectra of Probe Species in Complex Fluids 
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The objective of this paper is to correct a substantial, widespread error in parts of the quasielastic 
scattering literature. This error leads to entirely erroneous interpretations of quasielastic scattering 
spectra. The error, which is most prominent for interpreting spectra of dilute probe particles 
diffusing in complex fluids, arises from a valid calculation that is being invoked under circumstances 
in which its primary assumptions are incorrect. 

Quasielastic scattering from dilute probes yields the incoherent structure factor g (ls ^(q,t) — 
(exp(iqAx(t))} , with q being the magnitude of the scattering vector q and Ax(t) being the probe 
displacement parallel to q during a time interval t. The error is the claim that g^ ls '(q, £) for dilute 
probe particles uniformly reduces to exp(— q 2 ((Ax(t)) 2 )/2), regardless of the nature of the sur- 
rounding medium. If true, this claim would allow one to use quasielastic scattering to determine 
the time-dependent mean-square probe displacements in complex fluids. In reality, g^ l3 \q,t) is de- 
termined by all even moments ((Ax(t)) 2n ), n = 1, 2, 3, . . ., of the displacement distribution function 
P(Ax,t). Only in the very special case of monodisperse probes in a simple Newtonian solvent is 
g( l3 \q,t) entirely determined by {(Ax(t)) 2 ). Furthermore, the Langevin equation approach that 
ties g( ls \q,t) to {(Ax(t)) 2 ) also requires with equal certainty that g ( - la \q,t) relaxes as a simple 
exponential exp(— Ft), T being a time-independent constant. Contrariwise, if the spectrum is not 
a simple exponential in time, </ ls '(g, t) is not determined by ((Ax(t)) 2 ) . Several related subsidiary 
errors are discussed. 



I. INTRODUCTION 

Quasielastic scattering of light, x-rays, and neutrons 
has proven to be a powerful experimental tool for the 
study of complex fluids. For dilute macromolecule solu- 
tions, quasielastic scattering has extensive analytic appli- 
cations for particle sizing 11. With non-dilute solutions 
and more complex systems, quasielastic scattering reveals 
consequences of intermacromolecular and other interac- 
tions. 

In a substantial enhancement of the classical method, 
the fluid of interest is doped with trace concentrations of 
monodisperse probe particles. Probes have ranged from 
small molecules to micron-size colloidal particles. The 
diffusion of probe particles through the complex fluid is 
then observed. Early successful studies of probe diffusion 
in complex fluids using quasielastic light scattering were 
by Hallett and co-workersJH H|, who in 1974 and 1976 
observed the diffusion of polystyrene spheres through 
hyaluronic acid and dextran solutions, and compared 
probe diffusion with the rheological properties of their 
solutions. The subsequent four decades have seen an 
enormous extension of this approach^, including stud- 
ies of probes in highly viscous simple liquids [5|, poly- 
mer meltsH, chemically cross-linked gels p}, surfactant 
solutions [HT protein solutions and the interior of liv- 
ing cells [l0(. More recently, quasielastic x-ray scattering 
has been used to extend the range of distance scales over 
which diffusion can be observed [llj. 

Probe diffusion has also been studied by a series of 
other physical techniques, each technique being sensitive 
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to distinctive range of time and distance scales or other 
features of probe motion. For example, fluorescence cor- 
relation spectroscopy [T^], which by varying the probe 
concentration can measure both the self diffusion coef- 
ficient and the mutual diffusion coefficient of the labeled 
species [l3l flij. has in recent years been extensively used 
to measure tracer diffusion. Recent work using probe 
diffusion is sometimes termed microrheology, the term 
microrheology referring to a particular model [TBj for in- 
terpreting quasielastic scattering spectra. In some stud- 
ies, probe diffusion has been viewed as being of interest 
because it is a path to measuring the viscoelastic proper- 
ties of the solution. In other studies, probe diffusion has 
been viewed as being of interest because it measures so- 
lution properties that are not the same as the viscoelastic 
properties of the solution. 

A valuable complement to quasielastic scattering stud- 
ies of probe diffusion is provided by measurements on 
probes subject to external driving forces. The over- 
whelming bulk of the literature on driven probe mo- 
tion is provided by studies of capillary electrophoresis. 
The electrophoretic literature primarily emphasizes im- 
proving separations of different charged species. How- 
ever, electrophoretic separations often use solutions of 
neutral polymers as the support medium. At the same 
time as these experiments are performing separations, 
these experiments are also giving information on the 
dynamics of the neutral polymers 0, [llj]. A substan- 
tial literature exists on buoyancy- driven probe motion 
in the ultracentrifuge[I3, IS ■ In a few experiments, 
magnetic (l9l l2(ij or optical 21 1 tweezers were used to ex- 
amine oscillatory driven movements of probes. Tweezer 
experiments are particularly interesting because the ex- 
perimenter can separately control two of the three: drive 
force, drive frequency, and particle displacement. An al- 
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ternative complement to probe diffusion is provided by 
tracking probes in complex fluids in which the fluid itself 
is performing driven motion, e.g., shear (22j[. 

Quasielastic scattering of light and other photons 
is most commonly studied via correlation methods, in 
which the quantity measured directly is the intensity- 
intensity time correlation function 

S(q,t) = (I(q,T)I(q,T + t)), (1) 

Here q is the magnitude of the scattering vector, I(q, r) 
and I(q, r +t) are the scattering intensities over short 
time intervals near r and r + t, and the brackets (• • • ) 
represent an average. Scattering is said to be due to 
scatterers within the medium, scatterers generally being 
represented mathematically by points whose locations 
are known. Scattering from extended bodies, such as 
high-molccular-wcight polymer chains, is often treated 
by representing the scattering body as a series of scat- 
tering points whose relative positions are partly fixed. If 
the volume being observed is much larger than the vol- 
umes over which particle positions and displacements are 
correlated, quasielastic scattering corresponds to the in- 
termediate structure factor (or field correlation function) 
g (1 \q,t) via[H 

S(q,t) = A\gU(q,t)\ 2 + B. (2) 

In this equation A and B are constants determined by de- 
tails of the experimental apparatus; these constants have 
no effect on the time dependence. Homodync rather than 
heterodyne detection of the scattered light is assumed. 
The factorization of S(q,t) into g^(q,t) is sometimes 
termed the "Gaussian approximation" . This Gaussian 
approximation is not related to the Gaussian approxima- 
tion for the particle displacements as discussed below. 

The intermediate structure factor is in turn determined 
by the time-dependent positions of the scattering parti- 
cles via 

IN N \ 

9 (1) (q, *) = <££ cxp(*q • (r,(t + r) - r, (r))) ) (3) 

In this equation, sums on i and j proceed separately over 
all N particles in the system, while (t + r) and Tj (r) 
are the locations of scatterers i and j at times t + r and 
r, respectively. 

In applying eq [3j two particularly interesting experi- 
mental circumstances are described as mutual diffusion 
measurements and as probe diffusion measurements. In 
a measurement on a binary solvent: scatterer system, 
the scattering particles may be concentrated or dilute. 
Quasielastic scattering on such a system measures the 
mutual diffusion coefficient, which describes the diffusion 
of the scatterers down a concentration gradient [U [25|. 
Tracer diffusion experiments examine ternary solvent: 
matrix : probe systems. In these systems the matrix 
component may be dilute or concentrated, is substan- 
tially responsible for the system's rheological and other 



interesting properties, but is nearly optically inert. Con- 
versely, the probe (tracer) component is dilute, has vir- 
tually no effect on rheological and other properties of the 
solvent: matrix system, but dominates scattering by the 
ternary mixture. If matrix scattering is not entirely neg- 
ligible, there are established, reliable ways to isolate the 
probe scattering, based on spectral subtraction at the 
level of the field correlation function. 

Because probe particles very nearly do not interact 
with each other, the field correlation function for probe 
diffusion reduces (up to normalization constants) to the 
incoherent scattering function 

5 ( ls )(g,t) = (exp(zgAxW)). (4) 

with Ax(t) being the component parallel to q of Ar(i) = 
Yi{t + t) — r^(r). Probe motions perpendicular to q do 
not contribute to g^ ls ^(q, t). In moving from eq[3]to eq|H 
terms of eq[3]in which i ^ j were taken to average to zero, 
because the relative positions of dilute probes are very 
nearly uncorrelated. An expression formally identical to 
eq 2] describes diffusion measurements using pulsed- field- 
gradient nuclear magnetic resonance, though with this 
method q has an entirely different meaning, namely in 
the simplest case q = "/5g, where 7 is the gyromagnetic 
ratio, S is a pulse width, and g — dB/dz is the field 
gradient. 

The averages in eqs [3] and 0] may formally be phrased 
as averages over displacement distribution functions such 
as P(Ax,t), which gives the time-dependent probability 
that a scattering particle will dis plac e through Ax dur- 
ing time t. Two previous papers|26, |27( examined how 
g( ls \q,t) and g^~'{q,t) are actually related to the dis- 
placement distribution functions. The two prior papers 
were primarily concerned with establishing formal rela- 
tionships between dynamic structure factors and prob- 
abilities for scatterer displacement. The significance of 
these relationships for the interpretation of experimen- 
tal measurements was at most a secondary consideration. 
This paper focuses on interpreting experimental measure- 
ments. 

Section II of this paper presents the correct general re- 
lationship between g^ ls ^(q,t) and P(Ax,t). Section III 
discusses the special case of probe particles in a purely 
Newtonian fluid. Section IV notes experimental findings 
bearing on the relative significance of Sections II and III. 
Section V considers paths for interpreting probe diffusion 
spectra. Section VI treats the determination of P(Ax, t), 
relationships between g( ls \q,t) and trapping/hopping 
behavior, and, closing on a positive note, cases in which 
quasielastic scattering from diffusing probes, correctly in- 
terpreted, has given valuable information about complex 
fluids and the objects diffusing in them. 

II. GENERAL CASE 

This section summarizes what g^ ls \q,t) and g^(q,t) 
actually reveal about particle displacements. Extended 
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derivations have appeared previously in two earlier pa- 
pers, refs and [27], and are not repeated here. For 
probe diffusion, the intermediate structure factor is al- 
ways determined by the displacement distribution func- 
tion, namely the average in eq [4] can be written as 



d(Ax)e^q>(iqAx)P(Ax,t). (5) 



P(Ax, t) is taken to be properly normalized. 

On taking a Taylor series expansion of the exponential 
in powers of q, reflection symmetry, namely P(Ax,t) = 
P(—Ax,t), eliminates all terms odd in q. As a result, 
g( ls ^ (q, t) and its logarithm are necessarily power series in 
q 2 . The coefficients of the q 2n are generated by the even 
moments {(Sx) 2n ) of P(Ax,t). As shown previously [26j. 
the lead terms of an expansion for g( ls) {q,t) are 
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T x {t) is the random force, called random because in the 
Langevin model the values of T x (t) at different instants in 
time are uncorrelated. Within the model, T x cannot be 
predicted beyond stating that T x has certain statistical 
properties. 

The canonical literature treatment of the Langevin 
model as applied to quasielastic light scattering is the 
volume by Berne and Pecora |30|, notably their Section 
5.9. Berne and Pecora show that the Langevin model 
is appropriate for polystyrene spheres in water, on the 
time and distances scales observed by quasielastic light 
scattering. From the Langevin model and the require- 
ment that the system remains in thermal equilibrium, a 
series of conclusions about the statistical properties of 
the particle motion follow. In particular: 

(i) The mean-square average value of T x {fy must be con- 
sistent - the fluctuation-dissipation theorem - with 
the drag Coefficient / and the thermal energy fc^T. 



( 6 ) 

All even moments ((Sx) 2n ) are required for the complete 
expansion. 

It was early shown that quasielastic scattering from 
a binary solvent: macromolecule system determines 
the mutual diffusion coefficient, not the self diffusion 
coefficient [2J, |25[. Theoretical approaches to computing 
(q, t) and the mutual diffusion coefficient of non-dilute 
colloid solutions have historically followed routes very dif- 
ferent from the routes that are based on the displacement 
distribution function P(Ax,t). Only very recently [27l| 
was a solution for g^(q,t) in terms of displacement dis- 
tribution functions obtained. In this solution, the expan- 
sion of eq|3]was shown to require averages over two differ- 
ent displacement distribution functions, namely P(Ax, t) 
and a new distribution function Pi{Ax, t, R42). Pi is a 
two-particle conditional displacement distribution func- 
tion, in which Ax is the displacement of particle 1 during 
(0,t) given that the vector R42 from particle I to some 
particle 2 at time has a given value. 



III. SPECIAL CASE: PROBES IN SIMPLE 
NEWTONIAN LIQUIDS 

The earliest quasielastic scattering experiments were 
performed on dilute suspensions of monodisperse scat- 
tering particles in simple Newtonian solvents. Cummins, 
et al.'s results on polystyrene spheres in water [28| are the 
archetype. The resulting spectra were interpreted by in- 
voking a mechanical model for the motions of diffusing 
particles. The mechanical model was provided by the 
Langevin equation, which in one dimension is 



(li) The distribution P(Ax) of particle displacements 
Ace during a time interval At is the same for all 
time intervals (t, t + At). 

(iii) Velocity correlations are evanescent. For time steps 
appreciably longer than m/ f, which for Brownian 
particles is actually a quite short time, particle dis- 
placements in a series of time steps are very nearly 
independent from each other. 

Conclusion (ii) corresponds to the statement that x(t) 
is the sum of a series of identically-distributed random 
variables. Conclusion (iii) corresponds to the indepen- 
dent statement that the time evolution of x(t) is de- 
scribed by a Markoff process. In this very special case, 
the distribution of particle displacements is described by 
Doob's Theorem |29j|. Doob's theorem is closely related 
to the central limit theorem. Doob's theorem treats ran- 
dom processes such as Ax(t), while the central limit the- 
orem treats random variables. For the Langevin model, 
Doob's Theorem shows that the distribution of particle 
displacements is a Gaussian 

P(Ax) = (2Tr(Ax 2 )y 1/2 ex V (-(Ax(t)) 2 /2(Ax 2 )). (8) 

For this special case, the incoherent scattering function 
reduces to 



/ ls \q,t)=eM-q 2 ((^(t)) 2 )/2). 



(9) 



rn- 



d 2 x{t) 
' dt 2 



= ~f 



dx(t) 
dt 



(7) 



Here x(t) is a coordinate of the diffusing particle, m is 
the particle mass, / is the particle's drag coefficient, and 



Equation[9]is quite accurate for the systems considered in 
Ref. (3(| , namely highly dilute solutions of monodisperse 
objects in simple Newtonian solvents. 

However, Berne and Pecora [3(|, especially their Ap- 
pendix 5. A and Section 5.9 leading to their eq 5.9.6, also 
prove the other important consequence of the Langevin 
model and Doob's theorem, namely that the Langevin 
model determines the exact value of ((Ax(t)) 2 }. On the 
time and distance scales accessible to quasielastic scat- 
tering, the Langevin model requires 



((Ax(t)) 2 ) = 2Dt. 



(10) 
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Here ks is Boltzmann's constant and T is the absolute 
temperature. D = kgT/f is the diffusion constant, a 
quantity that does not depend on time. Time indepen- 
dence of D is required by the calculation, because D re- 
sults from a time integral over (0 < t < oo). 

Equations[9]and[T0]come as a package; they are equally 
consequences of the Langevin model. Correspondingly, 
Berne and Pecora show for diffusing monodispersc Brow- 
nian particles that the Langevin model requires that the 
field correlation function is a simple exponential 

g( ls \q : t) = exp(~q 2 Dt). (11) 

For unclear reasons - the literature error noted in the in- 
troduction - Berne and Pecora's entirely correct Chapter 
5 is being misread as proving that eq|H]is always correct, 
even when the time relaxation of g( ls \q,t) is not a sim- 
ple exponential. Berne and Pecora in fact prove exactly 
the opposite, namely the contrapositive of their result is 
that if the relaxation is not a single exponential, then 
the Langevin model must not be applicable to the sys- 
tem, and therefore invocation of the the Langevin Model 
prediction eq[9]is incorrect. 

Berne and Pecora's discussion refers purely and ex- 
clusively to the special case of particle motion that is 
described by the Langevin equation, in which case eqs 
I8KTT1 are all correct. This special case corresponds to 
most of the experiments that were of interest at the time 
that ref was written, namely applications of quasielas- 
tic scattering for particle sizing [31|. For diffusing dilute 
colloids, the t and q dependences of g^ ls ^(q,t) are pre- 
cisely as predicted by the Langevin model, in particu- 
lar g^ ls \q,t) ~ exp(— Ft) with T cx q 2 . The quasielas- 
tic scattering spectrum only leads to the time-dependent 
mean-square particle displacement if the spectrum is a 
pure exponential in q 2 t. In this special case, the mean- 
square displacement increases linearly with t, so that the 
short-time and long-time behaviors of g^ ls \q,t) are one 
and the same. 

If the decay of the field correlation function is not a 
simple exponential in tq 2 , then eq [8] and the Langevin 
model cannot not possibly describe how the scattering 
particles move. If eq [3] described the particle motions, 
then the spectrum would be a simple exponential. In 
systems in which the spectrum is more complex than a 
simple exponential, eq|H]is invalid. That is, if g^ ls \q,t) 
is not a simple exponential in t : log(t/ ls ) (q, t)) does not 
reveal the mean-square displacement of the particles. 

Why does eq |6] ever reduce to eq |9]? If and only if 
P(Ax, t) is a Gaussian in Ax, P(Ax, t) is entirely char- 
acterized by its second moment (Ax 2 ). For a Gaussian 
displacement distribution function, the higher moments 
of P(Ax, t) have values such that the coefficients of the 
higher-order terms (q 2n for n > 2) of cq [5] all vanish. 
For a Gaussian P(Ax, t), the only non-zero part of eq[6] 
is eq[9l This disappearance of the higher-order terms is 
unique to a Gaussian P(Ax,t). For any other P(Ax,t), 
the higher-order terms of eq |5] do not vanish. 



IV. EXPERIMENTAL FINDINGS 

What do experiments say about P(Ax, t) and 
g( ls \q,t)7 There are systems in which the Langevin 
model is adequate, namely dilute monodisperse particles 
suspended in simple Newtonian fluids. The Langevin 
model provides the solid foundation for particle sizing 
via quasielastic scattering [3l| . For probe diffusion in 
complex fluids, experiment provides a far more complex 
picture. Consider a few representative experiments that 
have determined PAx,t) or g^ ls \q,t) for probe particles 
in complex fluids. 

On relatively long time scales, P(Ax, t) is accessible via 
particle tracking methods. As examples, note experimen- 
tal studies by Apgar, et al.fHJ, Tseng and Wirtz(33l|. and 
Xu, et al.[34| on probe diffusion in glycerol, actin solu- 
tions and gels, and gliadin solutions. These authors used 
video recording and computer image analysis to track 
simultaneously the motions of large numbers of parti- 
cles. They report P(Ax, t) and (Ax(t) 2 ) in their systems. 
Probes in glycerol follow eqs 151 and [TUl For probes in vari- 
ous complex fluids, P(Ax, t) has decidedly non-Gaussian 
forms. In these systems, the mean-square displacement 
does not increase linearly in time. By direct measure- 
ment, eq [5] and [TU] are not uniformly correct for probes 
in polymer solutions. 

Quasielastic scattering spectra of probes in poly- 
mer solutions are often markedly non-exponential. For 
polystyrene latex sphere probes in hydroxypropylcellu- 
lose: water, this author and Lacroix[35| found stretched 
exponentials in time 

g {ls) (q,t) =aexp(-9t fl ) = aexp(-(t/r)' 3 ). (12) 

Here /3 is a scaling exponent while and r are prefac- 
tors. A series of papers by Streletzky and collaborators 
on the same chemical system (most recently, ref. l36l ) es- 
tablished by viewing a much wider range of delay times 
that g {ls) (q,t) is in fact a sum of two stretched exponen- 
tials in time. 

Finally, I note a very simple model system in which eqs 
IU and H] fail. The system is a dilute aqueous dispersion 
of polystyrene spheres, in which the spheres are of two 
different sizes. There are no sphere-sphere interactions. 
Each sphere individually performs Brownian motion as 
described by the Langevin equation. Therefore, for each 
sphere in the mixture, P(Ax, t) is a Gaussian in Ax and 
((Ax) 2 } increases linearly in time. For the mixture as a 
whole, the mean-square displacement averaged over all 
the particles must therefore also increase linearly with 
time, at a rate determined by a weighted average of the 
diffusion coefficients of the two sphere species. 

However, the mixture's field correlation function is a 
sum of two exponentials 

g^ (q, t) = A l cxp{-D iq 2 t) + A 2 exp{-D 2 q 2 t) (13) 

where Ai is the scattering amplitude of species i and Di 
is the diffusion coefficient of species i. Correspondingly, 
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the displacement distribution function for all the parti- 
cles in solution is a weighted sum of two Gaussians, one 
Gaussian for each sphere size. 

Suppose one were used eq|9]to determine ((Ax) 2 } of 
the sphere mixture. According to eq|H] 

((Ax) 2 ) = -\n(g^\q,t))/q 2 . (14) 

At short times, g^ ls \q,t) contains contributions from 
both species, and the {(Ax) 2 ) computed from eq[U]is the 
weighted average of the mean-square displacements of the 
two species. At long times, in eq[T3]the exponential corre- 
sponding to the more rapidly-moving species has decayed 
to zero, so the nominal mean-square displacement from 
eq 03] is determined by scattering from the larger, more 
slowly-moving spheres. For this simple bidisperse sphere 
suspension, eq HU asserts that the slope of {(Ax) 2 } de- 
pends on time, ((Ax) 2 ) increasing more rapidly at short 
times and more slowly at long times. The assertion is 
incorrect. In this system ((Ax) 2 ) increases linearly with 
t. At large times the smaller spheres continue to move, 
but they stop contributing to g^ ls \q, t) or to the nominal 
((Ax) 2 ) from eq[H 

Experiment thus demonstrates that neither eq[8]nor eq 
1111 is generally valid for probe diffusion in complex flu- 
ids. Even in a Newtonian fluid, a model system in which 
5 ( - ls ^(g,i) does not decay exponentially in time does not 
follow eqs [9] and [T4J a result that is exactly as required 
by Doob's theorem. Interpretations of quasielastic scat- 
tering spectra for probes in complex fluids, based on the 
Gaussian approximation of eq [HI are therefore incorrect. 
Interpretations of spectra of probes in complex fluids in 
terms of particle displacements are properly based on eq 
IH1 which correctly reflects the non- Gaussian displacement 
distribution function of real probes. 

V. INTERPRETATIONS OF QUASIELASTIC 
SCATTERING SPECTRA 

First, every single physical g^ ls ' viewed only as a func- 
tion of t corresponds to a system in which the mean- 
square particle displacement increases linearly with time. 
However, the correspondence is not unique. The same 
g^ ls \q, t) may also correspond to systems in which parti- 
cle thermal motions are more complex. In consequence, 
from a g^ 1 ^ (q, t) measured over a full range of times and 
a single q one cannot infer how the particle displacement 
depends on time. 

This result has a purely mathematical basis, namely 
that the field correlation function can always be repre- 
sented via a Laplace transform as 

/>00 

g {ls \q,t)= / A(T)exp(-Tt). (15) 
Jo 

Here T is a relaxation rate and A(r) is the contribu- 
tion of relaxations having decay rate T to g 1 - 1 ^ (q,t). So 
long as the system does not have relaxational modes that 



have negative amplitudes, A(T) is everywhere positive or 
zero. In this case, there is always a system having the 
same g( ls ^(q,t) as the system of interest, and in which 
((Ax(t)) 2 ) increases linearly in time. The system can 
be physically constructed as a solution of of polystyrene 
spheres of all different sizes. The composition of the mix- 
ture is determined by A(T): One adds to the mixture just 
enough polystyrene spheres having diffusion coefficient 
r/<7 2 so that their contribution to the scattering spectrum 
is ^4(r). For each sphere, {(Ax(t)) 2 ) increases linearly in 
time, so therefore ((Ax(t)) 2 ) of the mixture also increases 
linearly in time. Thus, an arbitrary physically-acceptable 
(i.e., A(T) > V r) form for A(T) corresponds as one 
non-unique possibility to a system in which ((Ax(i)) 2 ) 
increases linearly in time. 

It has repeatedly been found that g( ls ^(q,t) decays in 
time as the stretched exponential of eq [T2] If one inter- 
preted this time dependence by applying eq|9l one would 
conclude 

((Ax(t)) 2 H^. (16) 

In the common case /? < 1 , from eq [TH] one would infer 
that the mean-square particle displacement increases less 
rapidly at large times than would be the case for simple 
diffusion, a behavior that has been termed 'subdiffusive' 
motion. The inference is incorrect. A more reasonable 
interpretation for f3 < 1 is that diffusion in the complex 
fluid is created by modes having a range of relaxation 
times, some longer than others, the contribution of the 
slower modes to the spectrum becoming more important 
at longer times. 

It is not suggested that there cannot exist subdif- 
fusive motion. Such motion has unambiguously been 
observed experimentally. Amblard, et al.,[2l[ studied 
probe motion in f-actin solutions, using magnetic tweez- 
ers and video microscopy. Small-bead motion was diffu- 
sive; larger-bead diffusion was subdiffusive with f3 « 3/4. 
The non-classical drag forces for diffusive motion and for 
driven motion are the same, in the sense that the dis- 
placement under an applied force and the mean-square 
displacement for diffusion have the same linear or sublin- 
ear time dependences, depending on the probe size. 

There are sometimes suggestions that probe behavior 
should approach Stokes' Law and the Stokes Einstein- 
equation as the size of the probes is increased. Amblard, 
et al.'s experiments show precisely the opposite behavior. 
Deviations from simple diffusion are larger for the larger 
particles. Whenever particle motion is subdiffusive, the 
light scattering spectrum will not be a simple exponen- 
tial. The scattering spectrum will follow eq [51 showing 
that the relationship between g^ ls \q, t) and ((Ax(i)) 2 ) is 
a neither- way street. Just as one cannot in general cal- 
culate {(Ax(t)) 2 ) from g( ls \q,t), so also one cannot in 
general calculate g^ ls ^(q : t) from ((Ax(i)) 2 ), because all 
higher moments of P(Ax) are needed in order to calcu- 
late g( ls Hq,t). 
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VI. DISCUSSION 

The primary objective of this paper was to correct a 
widespread literature error, namely the false belief that 
that g' ls '(g,i) can in general be used to determine the 
mean-square displacement of probes in complex fluids. 
The belief appears to have arisen from a misreading of 
Berne and Pecora's excellent monograph, Chapter 5, in 
which Berne and Pecora discuss the motions of monodis- 
perse probe particles in simple Newtonian fluids by us- 
ing the Langevin model. The spectra of monodisperse 
Langevin-model particles are without exception single ex- 
ponentials. The calculation is correct, but does not refer 
to non-Newtonian fluids or polydisperse probe systems. 
Having corrected this error, I turn to several subsidiary 
points of misinterpretation. 

The functional form of P(Ax, t) can be inferred, at 
least approximately, from the angular dependence of 
</ ls ' (q,i), namely as seen in eq[S]the correlation func- 
tion g( ls \q,i) at fixed t is the spatial Fourier transform 
of P{Ax,t). If g( ls \q,t) is determined sufficiently ac- 
curately over an adequate range of q, an inverse spa- 
tial Fourier transform can take the experimenter back 
from g( la \q,t) to P(Ax,t). To the author's knowledge, 
this inversion has only been done for the ubiquitous 
polystyrene spheres in distilled water, for which g( ls ^ is a 
simple exponential exp(— Tt), while T is found to be accu- 
rately linear in q 2 , showing that P(Ax,t) in this system 
has the Gaussian form of eq[3 Measurements of the q- 
dependence of g^- 1 ^ for probes in complex fluids are less 
common, though note, e. g., Streletzky and Phillies (37|. 
These authors found spectra having multiple relaxational 
modes, some of which had relaxation rates that did not 
scale linearly in q 2 , proving that the modes did not corre- 
spond to Gaussian displacement distribution functions. 

Particle tracking is sometimes used to generate the 
simplified {Ax(t) 2 ), rather than the full P(Ax,t). The 
simplification is potentially hazardous, because if one has 
not determined P(Ax,t) one does not know if the parti- 
cle motion process corresponds to simple diffusion. Any 
physically reasonable P(Ax, t) has some second moment 
(Ax(t) 2 ), but if the form of P(Ax,t) is unknown, one 
cannot tell if it is meaningful to characterize P(Ax, t) 
by its second moment or by the corresponding nominal 
diffusion coefficient 



D(t) = (Ax(t) 2 )/2t 



(17) 



The notion that a diffusive process can be characterized 
by (Ax(t) 2 ) corresponds to particle diffusion in simple 
Newtonian liquids, in which D{t) as defined here is a 
constant independent of time. In a complex fluid, charac- 
terization of probe motion via measurement of the second 
moment (Ax(t) 2 ) must be expected to be inadequate. 

The assertion that the central limit theorem guaran- 
tees that P(Ax,t) is a Gaussian in Ax is sometimes de- 
scribes as the " Gaussian Approximation" . Experiments 
such as those summarized above prove that this asser- 
tion is incorrect. The central limit theorem (for random 



variables) and Doob's theorem (for random processes) 
are well-known. Where do their invocations go wrong? 
The central limit theorem and Doob's theorem are state- 
ments about the sum of a large number of identically 
distributed, independent random processes. As applied 
to the diffusion problem, the displacement Ax during ex- 
perimentally accessible times can be expressed as the sum 
of a large number of far smaller steps 8x, each taken dur- 
ing a far smaller time interval St. If a random variable 
Ax is the sum of a large number of identically distributed 
independent variables, i. e., if Sx(t) is a Markoff process, 
it must in general be the case that P(Ax, t) is a Gaussian 
in Ax. This rationale fails because it refers to a sum of 
independent random processes. The process that gener- 
ates 5x for probes in a viscoelastic fluid is not a Markoff 
process, because the system has memory. The central 
limit theorem and Doob's theorem are therefore not ap- 
plicable. For probes in complex fluids, the processes gen- 
erating the steps Sx are highly correlated, because the 
"random" forces that determine the Sx(t) are controlled 
by the shear modulus G(t), which in a complex fluid 
has a long correlation time. Correspondingly, the time 
correlation function of the random force (F X {Q)J- X (t)) is 
long-lived, not ~ S(t), and in the Langevin equation the 
friction force fx(t) is replaced with a memory function 
~fda(F x (0)r*(8))x(t-s). 

An alternative to the central limit theorem is the small- 
q approximation. The nominal idea in the small-q ap- 
proximation is that the rhs of eq|5]is a power series in q. 
If one went to sufficiently small q, one might hope that 
the q 2 term in the exponential would become dominant, 
so that eq[9]would approach being valid. This hope is not 
met. For the simplest case of a mixture of diffusing parti- 
cles, g^ ls ^(q, t) is in fact a power series in q 2 t. If one goes 
to smaller q, in order to keep fitting spectra equally ac- 
curately one needs to observe the same fractional degree 
of decay of g^ ls \q, t); one must therefore go out to longer 
times. At those longer times, the 0(q 4 ) terms are as sig- 
nificant as they were at larger q, smaller t, and the same 
q 2 t. Said differently, the coefficients of the correct Taylor 
series (in q) expansion of g^ ls \q,t) are time-dependent. 
In order for the lead term of the expansion to be domi- 
nant, the expansion must be limited not only to small q 
but also to to small t. If Ms large, no matter how small 
q has been made, the higher-order in q terms are as im- 
portant as the lower-order terms. Only at small q and 
small t is a single-term small-g expansion valid. The valid 
small-q expansion is 1 — q 2 (Ax(t) 2 ), which only describes 
the leading slope of g^ ls \q,t) at small times. 

Consider spectra described by eqU^l The exponential 
in eq [6] scales as q 2 or is a power series in q q , so there- 
fore 6 should also be a power series in q 2 , perhaps sim- 
ply by being linear in q 2 . Indeed, for probes in in some 
but not other water: hydroxypropylcellulose solutions, 



Streletzky [37] confirmed experimentally 9 ~ q over a 
wide range of q. If 9 were replaced with t _/3 , one would 
have 



-2//3 



(18) 
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f3 is often in the range 0.5-1.0, so r often depends on q 
as g~ 3±1 . If one interpreted r to be a relaxation time, 
the g-dependence from eq [18] would be strange indeed: 
The relaxation would occur more rapidly over large dis- 
tances (small q) than over short distances (large q). This 
strange g-dependence is simply an artifact of the choice 
of parameterization oi g( ls ^ (q, t) , and the identification of 
r as a relaxation time. In terms of eq |12[ {9, f3) provides a 
natural parameterization while (t, /3) is less transparent. 
If mean relaxation times are inferred from the spectral 
time moments 



t n g^ s \q,t)dt 



(19) 



of g^ 1 ^ (q, t), the choice of parameterizations in eq[T21has 
no consequences. The two paramaterizations of </ ls ) (q, t) 
must lead to the same (T n ). 

Spectra of diffusing probes showing two relaxations on 
very different time scales are sometimes interpreted in 
terms of caging and hopping relaxations. The notion is 
that the medium supplies regions of low potential en- 
ergy within which probes are free to move ("caging"). 
The regions are separated by barriers of high potential 
energy, across which probes only pass on rare occasion 
("hopping"). The short time-scale relaxation is said to 
correspond to caging, while the long time-scale relaxation 
is said to correspond to hopping. 

Computer simulation studies by Luo and Phillies and 
Luo, et al., test the caging-hopping interpretation [38l [39l j . 
These simulations represented Brownian particles mov- 
ing through a square lattice or a random glass of 
Lennard- Jones force centers. The force centers were im- 
mobile. Probe motions were generated via the Metropolis 
algorithm. These studies differed from some earlier work 
in that they determined not only time dependent mean- 
square displacements and effective diffusion coefficients 
but also obtained P(Ar, t) and g^ 1 ^ (q, t). By varying the 
nominal temperature, trapping, hopping, and hindered 
diffusion behaviors were obtained. At low temperatures, 
probe particles explored the volume of their traps; after 
a certain relaxation time (r 2 (t)) ceased to increase. At 
high temperatures, P(Ar, t) was nearly Gaussian, with 
(r 2 (t)) increasing linearly in time even at short times. 

Luo, et al., evaluated f/^ ls ^(g,i) for q^ 1 extending from 
a small fraction of the size of a single potential energy 
minimum out to distances substantially larger than a 
typical distance between force centers. At low and high 
temperatures, g^ ls '(q, t) showed nearly exponential relax- 
ations, though at small T and small q the relaxation fell 
to a non-zero baseline. The baseline was non-zero be- 
cause the particles were permanently trapped in small 
isolated volumes of the system. At intermediate temper- 
atures, relaxations were single-exponential at large q but 
double-exponential at small q. At the same intermediate 
temperatures, P(Ar, t) was radically non-Gaussian, with 
local maxima and minima created by local potential en- 
ergy minima, potential energy saddle points, and times 
required to traverse local energy maxima. 



Other, physically different, systems also give bimodal 
spectra. In contrast to Luo, et al.'s probes moving though 
a fixed matrix, in which relaxations are only bidisperse 
for some values of q, relaxations of dilute bidisperse sus- 
pensions are double-exponential at all q. An alternative 
model system in which monodisperse particles show sev- 
eral very different classes of relaxation behavior is shown 
by Glotzer, et al.'s[40] computer simulations of three- 
dimensional glasses, in which one finds distinct long-lived 
populations of slow and fast-moving particles, with the 
immobile particles in clumps and the rapidly moving par- 
ticles lying in thin ribbons. 

Thus, in order to distinguish between systems contain- 
ing species with two different dynamic behaviors, and 
systems in which there is local trapping with escapes from 
the traps at longer times, it is experimentally necessary 
to study g( ls \q,t) over a wide range of q. Observations 
at fixed q of double-exponential relaxations do not reveal 
whether one is seeing trapping with hopping, or whether 
the system is in some sense dynamically bidisperse. Fur- 
thermore, in the cases in which g^ ls \q,t) was observed 
by Luo, et al., to be very nearly the sum of two expo- 
nentials, P(Ar, t) on interesting distance scales had an 
elaborate dependence on r with multiple maxima and 
deep minima. The interpretation that a biexponential 
g( ls \q, t) must correspond to a P(Ar, t) that is a sum of 
two Gaussians, each with a mean-square width increas- 
ing linearly in time, is categorically disproved by Luo, et 
al.'s measured forms for P(Ar, t). 

Finally, the observation that quasielastic scattering 
does not determine the mean-square probe displacement 
certainly does not mean that probe diffusion is ineffective 
as an experimental technique. Probe diffusion measure- 
ments can certainly be used to obtain novel information 
about complex fluids. The richness of the revealed infor- 
mation corresponds to the depth with which models for 
probe motion are constructed. As a positive conclusion, 
two successful applications of probe diffusion are noted: 

(i) A long-time question in the study of surfactant so- 
lutions is the determination of the aggregation number 
n of surfactant molecules in micelles. One of many ap- 
proaches to this question has been to use quasielastic 
scattering to determine an effective hydrodynamic radius 
of the micelles. Perhaps after some hydrodynamic mod- 
eling to account for micelle shape, spherical micelles be- 
ing the simplest case, the measured diffusion coefficient 
can be transformed to an apparent hydrodynamic radius 
Th, to a hydrodynamic volume Vu, and (taking into ac- 
count the surfactant density and molecular weight) fi- 
nally to a nominal aggregation number. This procedure 
was criticized by Kratohvil 41], who noted that the hy- 



drodynamic volume of the micelle might well include sol- 
vent molecules rather than being composed of pure sur- 
factant. Probe diffusion experiments prove that Kratovil 
was correct. The diffusion of probe particles through 
micellar solutions is retarded by hydrodynamic and di- 
rect interactions between the micelles and the probe par- 
ticles. The degree of retardation is determined by the 
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volume fraction of micelles in the solution. By combin- 
ing quasielastic scattering measurements on surfactant 
solutions and on surfactant-probe mixtures, quasielastic 
scattering has been used to determine the size, volume 
fraction, and thus number density of micelles in solu- 
tion, leading to determinations of the micellar aggrega- 
tion number and, independently, the (substantial) degree 
of hydration of micelles, as seen in studies by Phillies and 
collaborators [42| . 

(ii) Diffusion of mesoscopic probe particles in polymer 
solutions is not Stokes-Einstcinian. D is not determined 
by the macroscopic viscosity 77. Therefore, one cannot 
use the Stokes-Einstein equation for sphere diffusion 



(where fcg is Boltzmann's constant, T is the absolute 
temperature, rj is the solution viscosity, and R is the 
sphere radius) to determine the size of probe particles in 
polymer solutions. However, by using probes of known 
size, Ullmann and Phillies [43j were able to quantitate the 
degree of failure of the Stokes-Einstein equation for their 
polymer solutions, allowing them to measure the size of 
unknown probe particles in the same solutions. This ap- 
proach permitted a quantitative study of the extent of 
polymer adsorption onto particles chosen for their abil- 
ity to bind polymers in solution. 
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